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Abstract: We propose a discrete model for how opinion about a given phe- 
nomenon, about which various groups of a population have different degrees of 
enthusiasm, such as fanaticism and extreme social and political positions, includ- 
ing terrorism, may spread. The model, in a certain limit, is the discrete analogue 
of a deterministic continuum model suggested by others. We carry out extensive 
computer simulation of the model by utilizing it on lattices with infinite- or short- 
range interactions, and on (directed) Barabasi-Albert scale-free networks. Several 
interesting features of the model are demonstrated, and comparison is made with 
the continuum model. 



1 Introduction 

Given the current political climate around the world, and the rise of extreme 
ideologies everywhere, from the Middle East to Africa and Western Europe, 
models that can provide insight into how such idealogies may spread in a 
society are clearly of great interest. In particular, given that, (1) the phe- 
nomenon of globalization has made interaction between people of various 
nations around the globe much easier than two decades ago, and (2) the fact 
that although extreme ideologies are usually advocated by very small fringe 
groups, but yet they continue to survive, it is important to understand the 
role of these two factors on the opinion about such antisocial behavior as 
terrorism. The goal of the present paper is to suggest a model to study this 
problem, and understand its implications. 

Some simple models for terrorism or extreme opinions appeared years 
ago in the physics PQ and sociological 2j literatures. The present work was 
motivated by an article on bioterrorism [3 , but the methods that we describe 
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and utilize can also apply to opinion dynamics regarding, for example, the 
latest "Star Wars" movie, fashion, a political candidate running for office, or 
other questions and opinions with varying degrees of enthusiasm. Thus, we 
do not even try to define "terrorism" here, as the model that we consider is 
generic. 

The population in our model consists of four parts, G, S, E, and F corre- 
sponding, respectively, to the general, susceptible, excited, and fanatic parts 
of the population. For simplicity, hereafter we use the same letters to denote 
the fractions of the total population belonging to each group. Members of the 
population can be convinced by acquaintances from the S, E, and F groups 
to move from the G group to S; from there by the E and F groups to change 
to E, and from there by F to change to F. Moreover, members of each of the 
three groups S, E, and F can change their status and go back directly to the 
G group. The dynamics of a model based on such a partitioning of a popu- 
lation has been treated in the continuum limit [3 by deterministic nonlinear 
differential equations, depending only on the total fractions G, S, E, and F. 
The continuum model can provide mathematically sufficient conditions for 
terrorism, or any other opinion about a certain subject, to die out at long 
times, implying that in the long-time limit everybody will belong to the G 
group , while the fractions S, E, and F shrink to zero which, when it comes 
to terrorism, is a good omen for the world. 

However, as is well-known in the statistical physics of complex systems, 
deterministic continuum models represent mean-field approximations to the 
actual problem which, although they allow for development of mathematical 
proofs for the existence or nonexistence of certain phenomena and provide 
us with a first guide, they are also unreliable. Such models cannot take into 
account the effect of fluctuations on the phenomena. In addition, such models 
cannot take into account the effect of the internet, fax machines, and satellite 
television which have made long-range interactions between people on very 
large scale possible. For example, in a phenomenon somewhat close to the 
present problem, deterministic differential equations predicted jl] extinction, 
whereas proper discrete simulations on a square lattice did not predict jl] 
the same phenomenon. 

Therefore, the goal of the present paper is to carry out extensive simu- 
lation of a discrete model opinion dynamics which, in a certain continuum 
limit, becomes similar to a deterministic model proposed by Ref. |3]. We uti- 
lize Monte Carlo simulations of a population of individuals. Such simulations 
may be called agent-based outside physics, but are used in physics since half 
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a century. 

The plan of the paper is as follows. We first describe the deterministic 
continuum model which is based on a set of nonlinear differential equations. 
We then describe the discrete model which is developed by putting individ- 
uals on a two-dimensional lattice, but the individuals can still be influenced 
by all other individuals. Then, we restrict the influence to nearest neigh- 
bours. Finally, we replace the regular 2D lattice by a scale-free network of 
Barabasi- Albert type jSj. The main point of the paper is not testing whether 
the model can provide quantitative predictions; rather, we deal only with the 
methods and how to implement them realistically. In particular we follow 
|3] in assuming that ifS = E = F = 0aX some moment, then these three 
quantities stay zero forever. Thus, one simulation corresponds to the opinion 
dynamics following one external event and does not include new external 
events to cause S, E, and F to become non-zero again. 

Differential equations 
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Figure 1: Differential equations: Fractions of susceptible, excited and fanatic 
opinions (from top to bottom) at 71 = 0.4: Fanatics die out. For 7! = 0.25, 
fanatics remain; for 7! = 0.55, excited and fanatic opinions die out, and for 
71 = 1.05, also the susceptibles die out. 
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2 The Deterministic Continuum Model 



The fractions G, S, E, and F in the population of agents having the cor- 
responding opinions, with C = S + E + F = 1 — G, change with time t 
as: 



Without loss of generality we set /?i = 1 since, otherwise, it can be absorbed 
in the time scale. We also set 72 = 7i- Nevertheless, we still have not only the 
four parameters 02, fa, 71, and 73, but also the three initial concentrations 
E(0), S(0), and F(0), which are relevant due to the nonlinearity of the 
continuum model. We use, in general, one million people and fa = 0.5, fa = 
0.5, and 73 = 0.20, starting with S = E = F = 0.2, which will be used 
throughout the paper. Then, we check for changes in the behaviour when we 
vary 71. 

Figure 1 illustrates the behaviour: As predicted in for 71 > fa (= 1 
here), only the general population remains; for decreasing 71 first, S also 
survives, then does E. Finally, for 71 = 0.25 F also survives, such that all 
the four groups, G, S, E, and F remain present in the population. 

3 Averaged Lattice 

Now we put the agents onto a 1000 x 1000 square lattice, half of whom is 
selected randomly to be fixed as empty. The fractions G, S, E, and F now 
refer to the filled half of the lattice, i.e., they are fractions of the population 
and not of the lattice size and, thus, still add up to unity. We try to follow 
closely the above set of equations, Eqs. (1), by the following rules for each 
time step t — > t + 1 (simultaneous updating of all agents): 

G becomes S with probability faC; 
F becomes G with probability 73; 

E becomes G with probability 72 and F with probability faF/C, and 
S becomes G with probability 71 and E with probability fa(E + F)/C. 

These changes are simulated by first looking at the decay through 7 and then 
at the radicalisation through (3. Therefore, it is possible that, e.g., an E first 
becomes G and immediately thereafter changes opinion to F. 



dS(t)/dt = faCG - p 2 S(E + F)/C - 7i5 
dE(t)/dt = [3 2 S(E + F)/C - (3 3 EF/C - l2 E 
dF(t)/dt = (3 3 EF/C - y 3 E . 



(la) 
(lb) 
(lc) 
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Figure 2: Averages on lattice: Fractions of susceptible, excited and fanatic 
agents (from top to bottom) at 71 = 0.4: Fanatics die out. For r y 1 = 0.2, 
fanatics remain; for 71 = 0.6, excited and fanatic agents die out, but even 
for 71 = 0.99, some susceptibles remain. 



Again we set f3\ — 1, 72 = 7i- Figure 2 looks similar to Fig. 1 which is 
not surprising since each agent is affected by all other agents, which is the 
limit that in statistical physics is described by a mean-field approximation. 
However, we always find for probability 71 < 1 some susceptible people, 
which is in contrast to their extinction by continuum model. Only in the 
unrealistic limit 71 = 1 do they die out. 

The reason for this persistence of susceptibles can be understood as fol- 
lows: For large enough 73, opinions E and F die out soon. Then, we have 
the differential equation for S = 1 — G as a simplification of Eq. (la): 

dS(t)/dt = - S)S - 7l S ~ (A - 71)5 (2) 

for small S, giving an exponential decay towards zero for 71 > (3\. For the 
Monte Carlo approach, G becomes S with probability (3\S and S becomes G 
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with probability 71. Equilibrium thus requires, for small S and thus G near 
unity, that: 



7l S = GfcS = &(! - S)S 



or 



S = lit Pi 



(3) 



which gives S — 1 — 71 for our choice /?i = 1. Only for 71 > /3i would S 
become zero, which is not possible if (3± — 1 since 71 is a probability for the 
Monte Carlo approach and no longer a rate which could also be larger than 
one. Putting back j3\ as a free parameter set equal to 0.5, everybody soon 
returns to the general population, S = E = F = 0, for 7x = 0.6. 



Nearest neighbours 
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Figure 3: Neighbour interactions on lattice: Fractions of susceptible, excited 
and fanatic agents (from top to bottom) at 71 = 0.4: excited and fanatic 
agents die out together. Both groups also die out at 71 = 0.9, while at 
71 = 0.1 the fanatics die out soon whereas excited agents decay very slowly. 
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4 Nearest Neighbour Interactions 



Now, we simulate a proper lattice population where only nearest neighbours 
influence each other. Thus at each time step every agent selects randomly 
one of the four nearest neighbours as a discussion partner and then follows 
these rules (again C = E + S + F): 

G becomes S with probability fti, if neighbour is S, E, or F; 

F becomes G with probability 73; 

E becomes G with probability 72, and 

it becomes F with probability fa/C, if neighbour is F; 
S becomes G with probability 71, and 

it becomes E with probability /fe/C, if neighbour is E or F. 

Thus, no agent is convinced by an empty neighbour to change opinion. 
Again we set (3i = 1, 72 = 7i- Since the behaviour of the population now 
depends on the single opinions and not only on their sum over all lattice sites, 
ordered sequential updating with helical boundary conditions was used. 

Figure 3 shows that, differently from Figs.l and 2, both E and F, and 
not only F, decay rapidly to zero; the susceptibles remain. 

5 Networks 

Human beings are not trees in an orchard, sitting on a square lattice and 
having k = 4 neighbours each. The previous section used a half-empty 
square lattice such that the number of neighbours varied between k = 
and k = 4. In reality, some people have many friends and some only few. 
Such social relationships are much better described by scale-free |H] networks 
of Barabasi-Albert type, where the number of people having k neighbours 
decays as 1/k 3 for not too small k. Moreover, real terrorism obeys a power 
law [7j. 

To grow such a network we start with four people all connected with each 
other and with themselves. Then, one after the other more sites are added, 
and each new site selects four people of the existing network as neighbours 
from whom to take advice. This selection is not random but proportional to 
the current number k of neighbours of that person: Powerful people attract 
more "friends" than powerless people. 

In this standard network, we then use directed opinion links. This means 
each person takes advice only from those four people whom the person se- 
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Scale-free network 
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Figure 4: Barabasi- Albert network: Now S, E, F all die out for 71 = 0.4, 0.9 
and (more slowly) 0.1. 



lected when joining the network; the same person ignores the opinions of 
those who joined later and selected this person as advisor. Directed networks 
have been used before for, e.g., the Ising models [S] and opinion dynamics 0. 
A computer program was listed in [TU]. Again, ordered sequential updating 
was used. 

Figure 4 shows that in the present model everyone becomes normal: S = 
E = F = after sufficiently long time, differently from the results obtained 
with the square lattice of the previous section. In contrast to the square 
lattices for 50 x 50 up to 20, 000 x 20, 000, larger networks needed a slightly 
longer time for E and F to decay; see Fig. 5. With different parameters, 
also quite complicated dynamics can be found; see Fig. 6. 

If we use symmetric instead of directed connections in these networks, now 
containing only 100,000 people, then for 71 = 0.1 all four groups survive; for 
71 = 0.4 the F die out; for 71 = 0.6 also the E and for 71 = 0.9 we end up 
with S = E = F = (not shown). Similar effects are seen if, initially, the 
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Figure 5: Barabasi-Albert network: Decay of excited fraction E becomes 
slower for network size increasing from 1000 (leftmost curve) to 30,000, then 
1 million and finally (rightmost curve) 30 million. We averaged over 10 
samples. 

first four people are fanatic and all others have opinion G. 

Thus, similarly to [10J for Ising magnets and differently from |9 for opin- 
ion dynamics, the directed network structure gives very different results com- 
pared to the undirected case. 



6 Conclusions 

In summary, the model of [3] depends somewhat on the various changes in 
the underlying connections which we had introduced. Overall we feel that we 
confirmed the main conclusions from differential equations 3J: Depending on 
parameters like 71, fanatics and/or excited agents survive or die out. 
DS thanks F. Bagheri-Tar for crucial help to survive in Los Angeles. 
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Figure 6: Barabasi-Albert network: Complicated dynamics at 7i = 72 = 
73 = 0.02 and initial S = E = F = 0.3. The last distribution remains stable 
for much longer times. 
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